Creating vortons and three-dimensional skyrmions from domain wall annihilation 
with stretched vortices in Bose-Einstein condensates 
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We study a mechanism to create a vorton or three-dimensional skyrmion in phase-separated two- 
component BECs with the order parameters ^i and ^2 of the two condensates. We consider a 
pair of a domain wall (brane) and an anti-domain wall (anti-brane) stretched by vortices (strings), 
where the ^2 component with a vortex winding is sandwiched by two domains of the ^1 component. 
The vortons appear when the domain wall pair annihilates. Experimentally, this can be realized 
by preparing the phase separation in the order ^1, ^2 and ^1 components, where the nodal plane 
of a dark soliton in ^1 component is filled with the ^2 component with vorticity. By selectively 
removing the filling ^2 component gradually with a resonant laser beam, the collision of the brane 
and anti-brane can be made, creating vortons. 

PACS numbers: 03.75.Lm, 03.75.Mn, 11.25.Uv, 67.85.Fg 



I. INTRODUCTION 

Quantized vortices are one of remarkable conse- 
quences of superconductivity and superfluidity. In multi- 
component superfluids and superconductors, there ap- 
pear many kinds of exotic vortices. When a vortex of 
one condensate traps another condensate inside its core, 
a supercurrent or superflow of the latter can exist along 
the vortex line. Such a vortex is called a superconducting 
or superflowing cosmic string in cosmology [1 . Because 
of the Meissner effect, superconducting strings exclude 
magnetic fields like superconductive wires, so that they 
are proposed to explain several cosmological phenomena 
related to galactic magnetic fields. When a supercon- 
ducting string is closed and "twisted", i.e., when the sec- 
ond condensate inside the string core has a non-trivial 
winding along the string loop, the supercurrent persis- 
tently flows along the loop and makes it stable. Such a 
twisted vortex loop is called a "vorton", a particle-like 
soliton made of a vortex |21 [3] . While vortons were dis- 
cussed in ^He superfluids [4l, they are considered to be 
a candidate of dark matter, and a possible source of ul- 
tra high energy cosmic ray. There have been a lot of 
study about their stability, interaction, and applications 
to cosmology [5 . 

On the other hand, three-dimensional (3D) skyrmions 
are topological solitons (textures) characterized by the 
third homotopy group tts (5/7(2)) :^ Z in a pion effec- 
tive field theory. Skyrmions were proposed to be baryons 
[6]. Since their proposal, the skyrmions have been stud- 
ied for a long time about their stability, interaction, and 
applications to nuclear physics [7 . 

Both 3D skyrmions and vortons have been fascinat- 
ing subjects in high energy physics and cosmology for 
decades, and a lot of works have been done already, but 



they have yet to be observed in nature. On the other 
hand, these topological excitations, 3D skyrmions J8HT4] 
and vortons [HI [16], can be realized in Bose-Einstein 
condensates (BECs) of ultracold atomic gasses. More- 
over 3D skyrmions and vortons have been shown to be 
topologically equivalent in two-component BECs [H [9] . 
BECs are extremely flexible systems for studying soli- 
tons (or topological defects) since optical techniques can 
be used to control and directly visualize the condensate 
wave functions [17]. Interest in various topological de- 
fects in BECs with multicomponent order parameters has 
been increasing; the structure, stability, and creation and 
detection schemes for monopoles [TSHU], knots [22] and 
non-Abelian vortices [2T have been discussed [24]. 

It is, however, still unsuccessful to create vortons and 
3D skyrmions experimentally, although the schemes to 
create and stabilize them have been theoretically pro- 
posed J8til4] . In the present study, we propose how 
to create vortons or 3D skyrmions in two-component 
BECs from domain walls and quantized vortices. Spe- 
cific examples of the system include a BEC mixture of 
two-species atoms such as ^'^Rb-'^^K [25^ or ^^Rb-^'^Rb 
[26] , where the miscibility and immiscibility can be con- 
trolled by tuning the atom-atom interaction via Fesh- 
bach resonances. Here, the domain wall is referred to as 
an interface boundary of phase-separated two-component 
BECs. Although the interface has a finite thickness, 
the wall is well-defined as the plane in which both com- 
ponents have the same amplitude. Since a description 
of two-component BECs can be mapped to the 0(3) 
nonlinear sigma model (NLaM) by introducing a pseu- 
dospin representation of the order parameter [27^^291, 
the resultant wall-vortex composite soliton corresponds 
to the Dirichlet(D)-brane soliton described in Refs. [30] - 
34 , which resembles a D-brane in string theory 



Such a D-brane soliton has been already numericahy con- 
structed by us in two- component BECs [38]. We have 
found that these composite solitons are energetically sta- 
ble in rotating, trapped BECs and are experimentally 
feasible with realistic parameters. Similar configuration 
has been also studied in spinor BECs [39] . 

A brane-antibrane annihilation was demonstrated to 
create some topological defects in superfluid ^He [40] . 
However, a physical explanation of the creation mech- 
anism of defects still remains unclear. The intriguing 
experiment that mimicked the brane-antibrane annihila- 
tion was performed in cold atom systems with the order 
parameters ^i and ^2 of two-component BECs by An- 
derson et at. [41] . They prepared the configuration of the 
phase separation in the order ^1, ^2 and ^1 components, 
where the nodal plane of a dark soliton in one component 
was filled with the other component. By selectively re- 
moving the filling component with a resonant laser beam, 
they made a planer dark-soliton in a single-component 
BEC. Then, the planer dark soliton in 3D system is dy- 
namically unstable for its transverse deformation (known 
as snake instability) [41 , which results in the decay of 
the dark soliton into vortex rings. In the two-component 
BECs, we have numerically simulated brane-anti-brane 
annihilations, which resulted in vortex loops [^. 

In this paper, we consider a junction of a D-brane soli- 
ton and its anti-soliton, namely a pair of a domain wall 
and an anti-domain wall stretched by vortices. We give 
an approximate analytic solution for a pair of the D- 
brane and anti-D-brane in the 0(3) NLaM. We show 
that this unstable configuration decays into a vorton or 
a 3D skyrmion, instead of an untwisted vortex ring [41] 
in the case without stretched vortices. Experimentally, 
this can be realized by preparing the phase separation 
in the order ^1, ^2 and ^1 components, and rotating 
the intermediate ^2 component. By selectively remov- 
ing the filling ^2 component gradually with a resonant 
laser beam, the collision of the D-brane and anti-D-brane 
can be made, to create vortons. 

This paper is organized as follows. In Sec. [TT] we 
present the Gross-Pitaevski energy functional of two- 
component BECs, and rewrite it in the form of NLaM. 
In Sec. |III[ after constructing a phase separation, i.e., a 
domain wall configuration in NLcrM, we consider a pair 
of a domain wall and an anti-domain wall. We discuss 
a creation of vortex in two dimensions, and a creation 
of vortex loops in three dimensions after a pair annihila- 



tion of the domain walls. In Sec. IV, we consider a pair 
of a domain wall and an anti-domain wall with vortices 
stretched between them. We show that when a vortex 
loop encloses n of the stretched vortices, the phase of the 
^2 component winds n times, i.e., it is a vorton with n 
twist. We also confirm a vorton with n = 1 is topologi- 
cally equivalent to a 3D skyrmion. Sec. [V] is devoted to 
a summary and discussion. 



II. SYSTEM 
A. Gross-Pitaevski energy functional 

The order parameter of two-component BECs is 

* = (*1,*2), (1) 

where 

% = v^e*^^ (j = l,2) (2) 

are the macroscopically occupied spatial wave function 
of the two components with the density pj and the phase 
Oj . The order parameter can be represented by the pseu- 
dospin 

s = (si , 52 , S3) = (sin cos (^, sin sin (/), cos 0) (3) 

with a polar angle = cos~^[{pi — P2)/ p] and an az- 
imuthal angle (j) = O2 — Oi as 



^ = y^e*2 (cos-e *2^sm-e*2 



(4) 



where p = pi -\- p2 and Q = Oi -\- O2 represent the local 
density and phase, respectively [27 . 

The solutions of the solitonic structure in two- 
component BECs are given by the extreme of the Gross- 
Pitaevski (GP) energy functional 
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Here, rrij is the mass of the jth component and pj is 
its chemical potential. The BECs are confined by the 
harmonic trap potential 



Vj = -mj{u;lx''+u;y+u;y). 



(6) 



The coefficients ^n, ^22, and gi2 represent the atom- 
atom interactions. They are expressed in terms of the s- 
wave scattering lengths an and a22 between atoms in the 
same component and ai2 between atoms in the different 
components as 
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with 771 ~^ = mj -\- m'^ . The GP model is given by 
the mean-field approximation for the many-body wave 
function and provides quantitatively good description of 
the static and dynamic properties of the dilute-gas BECs 



B. Mapping to the nonlinear sigma model 

To derive the generahzed NLaM for two-component 
BECs from the GP energy functional ([5|, we assume 
mi = 7Ti2 = Tn and V\ = V2 = V - By substituting the 
pseudospin representation Eq.Q of ^, we obtain [27] 
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where we have introduced the effective superflow velocity 
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Veff ^ 

and the coefficients 
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(12) 



The coefficient Ci can be interpreted as a longitudinal 
magnetic field that aligns the spin along the x^-axis; it 
was assumed to be zero in this study. The term with the 
coefficient C2 determines the spin-spin interaction asso- 
ciated with 53; it is antiferromagnetic for C2 > and 
ferromagnetic for C2 < [27 . Phase separation occurs 
for C2 < 0, which we are focusing on. Further simpli- 
fication can be achieved by assuming that V = and 
the total density is uniform through the relation p = p/g 
where 9 = 911 = 922 and p = pi = P2^ and that the 
kinetic energy associated with the superflow Veff is neg- 
ligible. Although the assumptions Veff = and p=const. 
become worse in the vicinities of vortex cores or domain 
walls, this simplification does not affect on later discus- 
sions about the vorton nucleations based on topology. 

By using the healing length ^ = %/ ^/2mgp as the 
length scale, the total energy can reduce to 
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where M is the effective mass for S3. This is a well-known 
massive NLaM for effective description of a Heisenberg 
ferromagnet with spin-orbit coupling. 
Introducing a stereographic coordinate 
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we can rewrite Eq. ([T3| as 

^3 1,:^ „.|2 I A/r2L,|2 
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Here, u = (00) corresponds to the south (north) pole 
of the 5*^ target space. 



III. WALL - ANTI-WALL ANNIHILATION 

A. Domain walls 

For a domain wall perpendicular to the x^-axis, u = 
u{x^)^ the total energy is bounded from below by the 
Bogomol'nyi-Prasad-Sommerfield (BPS) bound as [30l 



E = d^x 



\diu T Mu\^ ± M{u*diu + udm*) 



(l + H^)^ 



(17) 



by the topological charge that characterizes the wall: 



T. = Mi d3,!^:gi!i±^^ 



(1 + H2) 



(18) 



where di denotes the differentiation with respect to x*. 
Among all configurations with a fixed boundary condi- 
tion, i.e., with a fixed topological charge T^, the most 
stable configurations with the least energy saturate the 
inequality (17) and satisfy the BPS equation 



d^u T Mu = 0, 



(19) 



which is obtained by |...| = in Eq. (17). This equation 



immediately gives the analytic form of the wall configu- 
ration 



u^{x^) 



,TM{x^-xl)- 



(20) 



The function u^ represents the domain wall with wall po- 
sition Xq and phase ^0 associated with (si, 52); this phase 
00 yields the Nambu-Goldstone mode localized on the 
wall, as in Fig. [l] The sign ^ implies a domain wall and 
an anti-domain wall. The domain wall can be mapped to 
a path in the target space as shown in Fig. [2|a). 



B. Wall anti-Avall annihilation 

As described in SecjTj we note that the intriguing ex- 
periment that mimicked the brane-antibrane annihila- 
tion was performed by Anderson et al. [41 . They cre- 
ated the configuration shown in Fig. |3J where the nodal 
plane of a dark soliton in one component was filled with 
the other component. By selectively removing the fill- 
ing component with a resonant laser beam, they made 
a planer dark-soliton in a single-component BEC. The 
dark soliton corresponds to the coincident limit of the 
two kinks in Fig. [sFa). It is known that the planer dark 
soliton in the 3D system is dynamically unstable for its 
transverse deformation (known as snake instability) [4T] . 
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FIG. 1: (Color online) A single domain wall in two-component 
BECs. (a) The amplitude of ^i for a domain wall, (b) The 
pseudospin texture of the single domain wall perpendicular 
to the x^-axis in real space. The arrows denote points in 
the target space S^. The gradient and interaction energies 
are localized around the wall, which is shaded schematically. 
The arrows on the wall imply the phase 00 which the wall 
possesses. 
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FIG. 3: (Color online) A pair of domain wall and anti-domain 
wall in two-component BECs. (a) The amplitude of ^i for 
the wall and anti-wall configuration, (b) The pseudo-spin 
texture of the wall and ant i- wall configuration in real space. 
The arrows denote the pseudo-spin, ^i (^2) is filled outside 
(between) the walls, where the other component is zero. In 
the upper (lower) region outside the walls, the phase of ^1 is 
fixed to be zero (tt). 





FIG. 2: (Color online) The S^ target space where the north 
and south poles are denoted by and 0, respectively. 

(a) The path connecting the north and south poles represents 
the map from the path in the domain wall in Fig. [lib) along 
the x^-axis in real space from x^ -^ —00 to x^ -^ +00. The 
path in the S^ target space passes through one point on the 
equator, which is represented by "^" in Fig. [l[b) in this 
example. In general, the U{1) zero mode is localized on the 
wall. 

(b) The path in the target space S^ for a domain wall and 
an anti-domain wall. The path represents the map from the 
path along the x^-axis from x^ -^ —00 to x^ -^ +00 in real 
space in Fig. [3|b). 



which results in the decay of the dark soliton into vortex 
rings. 

In our context, this experiment demonstrated the 
wall-anti-wall collision and subsequent creation of cos- 
mic strings, where the snake instability may correspond 
to "tachyon condensation" in string theory [44]. The 
procedure that removes the filling component can de- 
crease the distance R between two domain walls and 
cause their collision. The tachyon condensation can leave 
lower dimensional topological defects after the annihila- 
tion of D-brane and anti-D-brane. In our case of the 



phase-separated two-component BECs, the annihilation 
of the 2-dimensional defects (domain walls) leaves 1- 
dimensional defects (quantized vortices). 

Let us discuss this in two dimensions in more detail. 
Here U{1) zero modes of the wall and the anti-wall are 
taken to be opposite as in Fig. [31(b). The configuration is 
mapped to a loop in the S'^ target space, see Fig. [2] (b). 
This configuration is unstable. It should end up with 
the vacuum with the up-spin 0. In the decaying process 
the loop is unwound from the south pole in the target 
space. To do this there are two topologically inequivalent 
ways, which are schematically shown in (a) and (b) in 
Fig. [4j In real space, at first, a bridge connecting two 
walls is created as in (c) and (d) in Fig. Hi Here, there 
exist two possibilities of the spin structure of the bridge, 
corresponding to two ways of the unwinding processes. 
Along the bridge in the x^-direction, the spin rotates (c) 
anti-clockwise or (d) clockwise on the equator of the 5^ 
target space. Let us label these two kinds of bridges by 
"I" and "t", respectively. 

In the next step, a 'passage' through the bridge is 
formed as in (e) and (f) in Fig. [4J where the ground 
state, i.e., the up-spin state, is filled between them. 
The phase of the filling ^1 component through the pas- 
sage is connected anti-clockwise or clockwise [Fig. ffl (g) 
and (h)] Let us again label these two kinds of passages by 
"I" and "t"? respectively. In either case, the two regions 
separated by the domain walls are connected through a 
passage created in the decay of domain walls. Once cre- 
ated, these passages grow to holes in order to reduce the 
domain wall energy. 

Several holes are created in the entire decaying process. 
Let us focus a pair of two neighboring holes. Then, one 
can find a ring of a domain wall between the holes as 
shown in Fig. [5] Here, since there exist two kinds of holes 
(t and ^), there exist four possibilities of the rings, (a) 
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FIG. 4: (Color online) Decaying processes of the wall- ant i- wall pair. (a,b) The loop in the pseudo-spin space is unwound in 
two ways. (c,d) A bridge is created between the wall and anti-wall. In this process there are two possibilities of spin structure 
along the bridge. (e,f) The upper and lower regions are connected with a 'passage' through the bridge being formed. (g,h) The 
^1 component is filled inside the passage, and the phase of ^i component is connected anti-clockwise or clockwise. 



tl, (b) It, (c) tt and (d) || in Fig.js] In all the cases, the 
^2 component is confined in the domain wall rings. The 
phase of ^1 component has a nontrivial winding outside 
the rings of types (a) and (b), whereas it does not have 
a winding outiside the domain wall rings of types (c) 
and (d). Consequently, the domain wall rings of types 
(c) and (d) can decay and end up with the ground state 
0. However, the decay of the rings of types (a) and (b) 
is topologically forbidden; they are nothing but coreless 
vortices. 

In the 0(3) NLaM, the domain wall rings of types 
(a) of (b) are the Anderson-Toulouse vortices [45 , or 
lumps in field theory [46,. The solutions can be written 
as {z = x'^ -\- ix^) 



^ = ^0 = y^ 



Z — Zj 



or u = uq 



(21) 



for a lump or an anti-lump, where Zi G G represent the 
position of the lump while and A^ G C* with |A^| and 
argAi representing the size and the U{1) orientation of 
the lump, respectively. In fact, one can show that these 
configurations have a nontrivial winding in the second 
homotopy group 7r2(6'^) o:^ Z which can be calculated 
from 
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The wall rings of (a) and (b) in Fig. p] belong to +1 and 



— 1 of 7r2(6'^), respectively. Namely they are a lump and 
an anti-lump, respectively. 

So far we have discussed two dimensional space in 
which domain wall is a line and a vortex is point-like. In 
three dimensions, domain walls have two spatial dimen- 
sions. When the decay of the domain wall pair occurs, 
there appear two-dimensional holes, which can be labeled 
by I or t in Fig. [OJ^a). Along the boundary of these two 
kinds of holes, there appear vortex lines, which in general 
making vortex loops, as in Fig. [OJ^b). This process has 
been numerically demonstrated ^2 . The vortex rings 
decay into the fundamental excitations in the end. 



IV. D-BRANE - ANTI-D-BRANE 
ANNIHILATION 



A. D-brane soliton 



The D-brane soliton by Gauntlett et. al. ^ can be 
reproduced in two-component BECs as follows [38]. For 
a fixed topological sector, vortices (a domain wall) par- 
allel (perpendicular) to the x-'^-axis, the total energy is 
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FIG. 6: (Color online) Decay of a domain wall pair in three 
dimensions. Domain walls are two-dimensional objects, (a) 
Two kinds of holes labeled by t and ^ are created after a 
domain- wall pair decay, (b) Along the boundary of these two 
kinds of holes, there appear vortex lines (denoted by dotted 
lines), which in general making vortex loops. 



FIG. 5: (Color online) Stable and unstable wall rings. (a,b) 
Stable wall rings. The domain wall phase winds once (the 
winding number is ±1) along the rings. Total configurations 
are 2d skyrmions with a non-trivial element ±1 of the second 
homotopy group 7r2. The phase ^i winds ±1 so that it is a 
vortex, while ^2 is filled inside the core. (c,d) Unstable wall 
rings. The domain wall phase does not wind (the winding 
number is 0) along the rings. The phase ^1 does not wind so 
that it cannot be a vortex, while ^2 is filled inside the ring. 
They decay into ground state (up pseudo-spin). 



where l32l 



bounded from below by the EPS bound as [3Ql [32l [33 l [38] 
E = I d-'x- 
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by the topological charges that characterize the wall and 
vortices: 
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Then, the least energy configurations with fixed topo- 
logical charges (a wall with a fixed number of vortices) 
satisfy the BPS equations 

diu tMu = 0, (S2 T id3)u = 0. (26) 

The analytic form of the wall-vortex composite solitons 
can be found {z = x'^ -\- ix^) 



u{x^,z) = u^{x^)u^{z), 



(27) 
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The function u^^ represents the domain wall with wall po- 
sition Xq and phase (pQ. The function Uy gives the vortex 
configuration, being written by arbitrary analytic func- 
tions of z; the numerator represents Ny.^ vortices in one 
domain (^1 component) and the denominator represents 
Ny^ vortices in the other domain (^2 component). The 
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v(2) 



positions of the vortices are denoted by z- and z- . The 
total energy does not depend on the form of the solution, 
but only on the topological charges as T^ = ±M or 
(per unit area), and T^ = 27rA^v (per unit length), where 
Ny is the number of vortices passing through a certain 
x^ = const plane. 

Figure [7| shows a D-brane soliton with the simplest 
wall-vortex configuration of Eq. (27). A vortex exists in 



X < and forms a texture, where the spin points down 
at the center and rotates continuously from down to up 
as it moves radially outward. The edge of vortex attaches 
to the wall, causing it to bend logarithmically as x^ = 
log \z\/M [Fig.lTFb) and (d)]. We can construct solutions 
in which an arbitrary number of vortices are connected to 
the domain wall by multiplying by the additional factors 
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zy I see Eq. (29)]; Fig. 8 shows a solution in which 



both components nave one vortex connected to the wall. 
In the NLcrM, the energy is independent of the vortex 
positions zj on the domain wall; in other words, there 
is no static interaction between vortices. 
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FIG. 7: (Color online) The typical D-brane soliton in two- 
component BECs. (a) Schematic illustration of the wall- 
vortex soliton configuration viewed from the length scale 
larger than the domain- wall width and the vortex core size. 
The two-component BECs ^i (x^ < 0) and ^2 [x^ > 0) 
are separated by the domain wall. A single vortex located 
at x^ < (^1 component) is connected to the domain wall. 
(b) The isosurface of S3 = for the solution Eq. ( 27 ) of the 
NLcrM, where M = 1, xj = 0, 00 = 0, A^^i = 1, A^^2 = 0, and 
z[ ■' =0. The corresponding spin textures s in the 2; = plane 
and y = plane are shown in (c) and (d), respectively. The 
magnitude of S3 is denoted by color. We have s = (0,0, —1) 
along the vortex core. 



B. Brane-anti-brane annihilation with a string 

We are ready to study a pair of a domain wall and an 
anti-domain wall stretched by vortices. An approximate 
analytic solution of the domain wall pair stretched by one 
vortex, which is schematically shown in Fig. loFa), can be 
given in the 0(3) NLaM as 

u{x^,z) = u^{x^)u^{z), (30) 

u^{z) = 1/z. (32) 

Here, xl and x^ {xl < x^) represent the positions of the 
wall and anti-wall, respectively, while 0i and 02 denote 
the phase of the wall and anti-wall, respectively. This 
solution is good when the distance \x\ — xW between the 
walls is large compared with the mass scale M~^ . For 
our purpose, the phases are taken as 0i = 02 + tt, which 
means that the ^1 component has a dark soliton when 
the intermediate ^2 component vanishes. The isosurface 
of S3 = and the pseudospin structure of this configu- 



FIG. 8: (Color online) The D-brane soliton to which two vor- 
tices attach, (a) Schematic illustration of the configuration in 
which each component has a single vortex connected to the 
wall, (b) The isosurface of S3 = for the solution Eq. (27) of 
the NLcrM, where M = 1, xj = 0, 0o = 0, iV^i = 1, A^^2 = 1, 
z\ — —2, and z\^ —2. The corresponding spin textures s 
in the 2; = plane and 2/ = plane are shown in (c) and (d), 
respectively. The magnitude of S3 is denoted by color. The 
wall becomes asymptotically flat due to the balance between 
the tensions of the attached vortices. 



ration are plotted in Fig. loFb) and (c), respectively. In 
order to avoid the logarithmic bending of the walls, one 
can use u^{z) in Eq^_j29| with N^^ — Ny^ instead of 
Eq. (|32|, as in Fig. [lO|The solution in Eqs. ([30|-([32| 
of the 0(3) NLcrM has a singularity at the midpoint of 
the vortex stretching the domain walls, as in Fig. [9J It 
is, however, merely an artifact in the NLcrM approxima- 
tion of p = const.; the singularity does not exist in the 
original theory without such the approximation, because 
p varies and merely vanishes at that point. 

Now let us discuss the dynamics of the wall-ant i- wall 
configuration. As in the case without a stretched string, 
the configuration itself is unstable to decay, and vortex- 
loops are created in the ^1 component. Since the ^2 
component is localized along the vortex core, the 53 = 
surface forms a torus (ring), where the region of 53 > 
is outside the torus whereas the region 53 < is inside 
it. As the phase of ^2 component inside the ring is con- 
cerned, the vortex loops are classified into 1) the un- 
twisted case [see Fig. pTj^a)] and 2) the twisted case [see 

Fig.jnjb)]. 

1) If the closed vortex- loop encloses no stretched vor- 
tices as the loop A in Fig. 10, the vortex-loop is not 



twisted, as in Fig. 11 ^a). Equivalently, the phase of the 
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FIG. 9: (Color online) (a) A pair of a D-brane (domain wall) and an anti-D-brane (anti-domain wall) stretched by a string 
(vortex) in two-component BECs. The branes are perpendicular to the x^-axis and the string is placed along the x^-axis. The 
arrows denote pseudo-spins. The ^i (^2) component is filled outside (between) the branes, where the other component is zero. 
In the upper (lower) region outside the branes, the phase of ^1 is fixed to be zero (tt). In the middle region, the phase of ^2 has 
winding around the vortex placed at the x^-axis. Accordingly, the pseudo-spin rotates once ( ant i-) clockwise at the endpoint of 
string on the (anti-)brane. The profile of the ^1 component along a line parallel to the x^-axis at (x^,x'^) 7^ represents two 
kinks in the left pannel, while the profile of the ^1 component along the x^-axis shows the coincident two kinks, i.e., a dark 
soliton. The dot in the center denotes the point (^1, ^2) = 0, which corresponds to a singularity in the NLcrM approximation 
(p = const.) in (b). (b) The isosurface of S3 = of an approximate solution in Eq. (30) with a domain wall and an anti-domain 
wall stretched by a vortex in the 0(3) NLcrM, where M = 1, xi = —3, x\ = 3, 01 = 0, 02 = tt. (c) The pseudo-spin texture of 
an approximate solution in Eq. (|3Q|). 




FIG. 10: (Color online) Loops in the wall- vortex systems. 
While the loop A yields an untwisted vortex ring in Fig. 11 'a), 
the loop B (C) yields a vorton, i.e., a vortex ring twisted once 
(twice). A vorton with twisted once is shown in Fig. 11 ^b). 



^2 component inside the ring is not wound. 

2) However, if the vortex-loop encloses n stretched vor- 
tices as the loops B and C in Fig. [lO| the vortex-loop is 
twisted n times. It implies that the phase of the ^2 com- 



ponent inside the ring is wound n times. A vortex-loop 
twisted once, which is nothing but a vorton with the min- 
imum twist, is shown in Fig. flTlb). The vertical section 
of the torus by the x^-x^ plane is a pair of a skyrmion 
(coreless vortex) and an anti-skyrmion (coreless vortex). 
Moreover, the presence of the stretched vortex implies 
that the phase winds anti-clockwise along the loops, as 
can be seen by the arrows on the top and the bottom 
of the torus in Fig. [TTj^b). When the 2D skyrmion pair 
rotate along the x^-axis their phases are twisted and con- 
nected to each other at the tt rotation. Note that the ze- 
ros of ^1 and ^2 make a link. Along the zeros of ^1 (^2)5 
the phase of ^2 (^1) winds once. The configuration is 
nothing but a vorton. 

It may be interesting to point out that this spin texture 
is equivalent to the one of a knot soliton [22^ .28^ 29j . 
i. e. a topologically nontrivial texture with a Hopf charge 
7r3(6'^) !:^ Z in an 0(3) NLaM. Mathematically, this fact 
implies that a vorton is Hopf fibered over a knot. 

Finally, to confirm a vorton creation of a domain wall 
pair annihilations, we show a numerical simulation of the 
time-dependent GP equation ifidt^ j = 5E /d^"", for the 
domain wall pair with a stretched vortex in Fig. [12] The 
numerical scheme to solve the GP equation is a Crank- 
Nicholson method with the Neumann boundary condi- 
tion in a cubic box without external potentials. The box 
size is 52.1^ x 52.1^ x 52.1^ with ^ = %/^miii. We pre- 




(a) 




FIG. 11: (Color online) (a) The pseudo-spin texture of an untwisted vortex-ring, and (b,c) the pseudo-spin and phase textures 
of a twisted vortex ring, i.e., a vorton, after the brane-anti-brane annihilation, (b) The pseudo-spin structure of the vorton. 
The torus divides the regions of ^i and ^2 which repel each other: ^1 (^2) are filled outside (inside) the torus. The vertical 
section of the torus by the x^-x^ plane is a pair of 2D skyrmions and anti-skyrmions. While they rotate along the x^-axis their 



pseudo-spins are twisted. This spin texture is equivalent to that of a knot [221 128| [29] . (c) The phase structure of the vorton. 
The arrows denote the phase of ^1 and ^2- The circle denotes the core of vorton where ^2 is filled and ^1 is zero. The phase 
of ^2 winds once along that circle. The square of the dotted line denotes a loop where ^2 is zero. The phase of ^1 winds once 
along that loop. Note that the loops of ^1 and ^2 are zero respectively make a link. Along the zeros of ^1 (^2), the phase of 
^2 (^1) winds once. 



pare a pair of a domain wall and an anti-domain wall at 
coincident limit with a vortex winding in the ^2 compo- 
nent. Here, for simplicity we put a cylindrically symmet- 
ric perturbation, which is expected to be induced from 
varicose modes of the string. Several holes grow after be- 
ing created, and there appear vortex loops. Although the 
holes appear asymmetrically because of the cubic bound- 
ary, the boundary effect is small in the center region and 
the initial perturbation causes a vortex loop there. The 
vortex loop enclosing the ^2 winding, which is nothing 
but a vorton, is created in the center of Fig. 12 ^c). 



C. 



Equivalence of the vorton to the 
three-dimensional skyrmion 



It has been already shown in [HI |9] that 3D skyrmions 
are topologically equivalent to vortons in two-component 
BECs. In this subsection, we show it in our context of 
the brane-anti-brane annihilations. 

In Fig. [13J the arrows denote the phase of ^1 along a 
large loop (of the square of the dotted line) going to the 
boundary where ^2 is zero, making a link with the vorton 
core. The left panel of Fig. [13] represents the phases 
of ^1 and ^2 of the vorton from the brane-anti-brane 
annihilation [see also Fig. ITTlb)]. This is topologically 
equivalent to the right panel of Fig. [13] Here we show 



that the phase structure of the right panel is that of a 
3D skyrmion. They are topologically isomorphic to each 
other. 

First, let us introduce the matrix U as 



^1 
^2 



^1 -^5 
^2 ^t 



U 



(33) 



with 



U-: 



*1 

*2 *! 



n 



being an element of an 5/7(2) group, when 



detU = 1*1 1 



1*2 



1. 



(34) 



(35) 



The GP energy functional given in Eq. ^ is not SU{2) 
symmetric in general. When the relations 



gil =922 =912, Ml =/i2 



(36) 



hold, the GP energy functional is SU(2) symmetric, and 
Eq. ( [35| ) holds (up to overall constant) [47 . 

Even when the GP energy functional given in Eq. (5| 
is not SU(2) symmetric, we can approximately consider 
the parametrization by U in Eq. (34). A rotationally 



symmetric configuration of a 3D skyrmion can be given 
byi 



U = exp i 



/«^- 



(37) 



with a function /(r) with the boundary condition 

/(r = 0) = mr, f{r = R) = 0, (38) 

where R is the size of the system. Here, n G Z is an 
element of the third homotopy group 7rs{SU{2)) :^ Z. 
In the polar coordinates (r, 6>,0), 



= (sin 6 cos ^, sin 6 sin 0, cos 6) . 



(39) 
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(a) 



(b) 



(c) 



FIG. 12: (Color online) A numerical simulation of a vorton creation. Surfaces are defined by ni — n2 = 0.18(/ii/^) (/ii = g = 
l^Ti — m — 1), while color represents the phase of the ^2 component, (a) First, we prepare a pair of a domain wall and an 
anti-domain wall at coincident limit with a vortex winding in the ^2 component, (b) Holes are created in the wall annihilation. 
(c) A vorton is created in the center. 
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FIG. 13: (Color online) The equivalence between the vorton 
and the 3D skyrmion. The arrows denote the phase of ^1 
along a large loop (of the square of the dotted line) going to 
the boundary where ^2 is zero, making a link with the vorton 
core. The left panel represents the configuration of the vorton 
from the brane-anti-brane annihilation [see Fig. |ll| -(a)], while 
the right panel represents the configuration of a 3D skyrmion. 
They are topologically isomorphic to each other. 



By using the formula, exp(z6n-cr) = cos 6+m-cr sin 6 for 
n^ = 1, the 3D skyrmion in Eq. (33) with U in Eq. (37) 
can be obtained as 



Along the x^-axis {0 = 0,7r), Eq. (40) becomes 

at l9 = 0, 
at ^ = TT, 



^1 
^2 



exp[-z/(r) 




^i\ _ /exp[z/(r)] 
^2 ) \ 



(42) 



Eqs. (41) and (42) show, in the case of n = 1, the phase 



structure of ^1 in the right panel of Fig. 13 



Second, let us study the phase structure of ^2- We 
consider the ring defined by ^ = 7r/2 and the radius r 
such that fir) = 7r/2. Along this ring, the ^i are 



*2 







(43) 



The ^2 component winds once along this ring, as in 



Fig. 11 ^c). This winding of the ^2 component originates 
from the winding in the brane-anti-brane configuration 
in Eq. (|32|. 



We thus have seen that the 3D skyrmion in Eq. ( 33 ) is 
topologically equivalent to a vorton. 



^1 
^2 



cos /(r) — i sin /(r) cos 
sin/(r) sin^e"*"^ 



(40) 



First, let us study the phase structure of ^1 of the 3D 
skyrmion in Eq. (40). At the boundaries and the origin, 
Eq. (40) becomes 



^1 
^2 

^2 



1 


(-ir 




at r = i?, 



at r = 0. 



(41) 



V. SUMMARY AND DISCUSSION 

We have studied a mechanism to create a vorton 
or three dimensional skyrmion in phase separated two- 
component BECs. We consider a pair of a domain wall 
and an anti-domain wall with vortices stretched between 
them. The ^2 component is sandwiched by the regions of 
the ^1 component, where the phase difference of ^I's in 
the two separated regions is taken to be tt. When the do- 
main wall pair decays, there appear vortex loops of the 
^1 component with the ^2 component trapped inside 
their cores. If a ^1 vortex loop encloses one stretched 



11 



vortex, it becomes a vorton. More generally, if the vor- 
tex loop encloses n of the stretched vortices, it becomes a 
vortex ring with the phase of ^2 twisted n times. We also 
have confirmed that the vorton (n = 1) is topologically 
equivalent to a 3D skyrmion. 

Experimentally this can be realized by preparing the 
phase separation in the order ^1, ^2 and ^1 components, 
and rotating the intermediate ^2 component. By selec- 
tively removing the filling ^2 component with a resonant 
laser beam, the collision of the brane and anti-brane can 
be made, to create vortons. 

Once created in the laboratory, one can study the sta- 
bility and dynamics of a vorton experimentally. The vor- 
ton will propagate along the direction perpendicular to 
the initial configuration of the branes. Therefore, to in- 
vestigate the dynamics of a vorton, we need to prepare 
a large size of cloud in that direction. In the case of 
untwisted vortex loop (usual loop, not a vorton), it will 
easily shrink and eventually decay into phonons if the 
thermal dissipation works enough. However, the vorton 
should be stable against the shrinkage and will propagate 
to reach the surface of the atomic cloud. Such a difference 
must be a benchmark to detect vortons in experiments. 

On the other hand, the thermal and quantum fiuctu- 
ations may make the vorton unstable. Our numerical 
calculations rely on the mean field GP theory. The topo- 
logical charge of a vorton is the winding of the phase of 
^2 along the closed loop (which is proportional to the 
superfiow along the closed loop). Since this topologi- 
cal charge is defined only in the vicinity of the vorton, 
there is a possibility that it can be unwound, once quan- 
tum/thermal fiuctuation is taken into account beyond 
the mean field theory. Quantum mechanically, such a 
decay is caused by an instanton effect (quantum tunnel- 
ing). This process also resembles the phase slip of su- 
perfluid rings. The vorton decay by the quantum and 
thermal tunneling is considered to be an important pro- 
cess in high energy physics and cosmology, since it will 
radiate high energy particles such as photons, which may 
explain some high energy astrophysical phenomena ob- 
served in our Universe. Therefore it would be important 
that one realizes vortons in laboratory by using ultra-cold 
atomic gases; it may simulate a vorton decay with emit- 
ting phonons quantum mechanically, beyond the mean 
field approximation. 

In this paper, we have mainly studied topological as- 



pects of the vorton creation using the NLaM approxima- 
tion. In order to study dynamics of topological defects 
beyond this approximation, we need a precise form of the 
interaction between the defects. An analytic form of the 
interaction between vortices was derived in the case of 
miscible (c2 > 0) two-component BECs, and it was ap- 
plied to the analysis of vortex lattices [48 . Extension to 
the case of the immiscible case (c2 < 0) focused in this 
paper will be useful to study the interaction between vor- 
tices attached to domain walls, and that between vortons 
and/or walls 

In our previous paper [38], we discussed that the do- 
main wall in two-component BECs can be regarded as a 
D2-brane, as the D-brane soliton J3Qti34] in field theory, 
where "Dp-brane" implies a D-brane with p space dimen- 
sions. This is because the string endpoints are electrically 
charged under U{1) gauge field of the Dirac-Born-Infeld 
(DBI) action for a D-brane [49 . In our context, the U{1) 
gauge field is obtained by a duality transformation from 
U{1) Nambu-Goldstone mode of the domain wall. Since 
the D-brane soliton [38 in two-component BECs, pre- 
cisely coincides with a Blon [50 , i.e., a soliton solution 
of the DBI action of a D-brane, the domain wall can be 
regarded as a D2-brane. 

On the other hand, it is known in string theory 04] 
that when a Dp-brane and an anti-Dp-brane annihilate 
on collision, there appear D(p — 2) branes. If we want to 
regard our domain wall as a D2-brane, the pair annihi- 
lation of a D2-brane and anti-D2-brane should result in 
the creation of DO-branes. Therefore, a discussion along 
this line leads us to suggest a possible interpretation of 
3D skyrmions as DO-branes, which are point-like objects. 
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